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Abstract 

We discuss the phase diagram of Wilson fermions in the vriQ-g 2 plane for two- 
flavor QCD. We argue that, as originally suggested by Aoki, there is a phase 
in which flavor and parity are spontaneously broken. Recent numerical results 
on the spectrum of the overlap Hamiltonian have been interpreted as evidence 
against Aoki's conjecture. We show that they are in fact consistent with the pres- 
ence of a flavor-parity broken "Aoki phase" . We also show how, as the continuum 
limit is approached, one can study the lattice theory using the continuum chi- 
ral Lagrangian supplemented by additional terms proportional to powers of the 
lattice spacing. We find that there are two possible phase structures at non-zero 
lattice spacing: (1) there is an Aoki phase of width Amo ~ a? with two massless 
Goldstone pions; (2) there is no symmetry breaking, and all three pions have an 
equal non-vanishing mass of order a. Present numerical evidence suggests that 
the former option is realized for Wilson fermions. Our analysis then predicts the 
form of the pion masses and the flavor-parity breaking condensate within the 
Aoki phase. Our analysis also applies for non-perturbatively improved Wilson 
fermions. 
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I. INTRODUCTION 



This paper concerns the phase diagram of Wilson fermions (and improved versions 
thereof) at non-zero lattice spacing, and the restoration of chiral symmetry in the con- 
tinuum limit. Some time ago, Aoki proposed a phase diagram in which there were regions of 
spontaneous flavor and parity breaking (which we refer to as Aoki phases) [|IJ. This sugges- 
tion, sketched in Fig. [1], provides a dynamical explanation for the masslessness of the pion on 
the lattice. Analytical and numerical support for this proposal has been given in Refs. []]-§[]. 

The validity of this picture has, however, been challenged recently by Bitar, Heller and 
Narayanan M. These authors argue that flavor and parity violation do not occur at non- 
zero lattice spacing, and that, by analogy with the massless continuum theory, one should 
interpret the Aoki phases as containing massless quarks and exhibiting spontaneous chiral 
symmetry breaking. Numerical results supporting this claim are given in Refs. |7|,|8|]. 

In this paper we present two arguments supporting Aoki's suggestion of spontaneous 
flavor and parity breaking at non-zero lattice spacing. Our first observation is that the 
numerical results of Refs. ]/]J§ are not only consistent with Aoki's picture, but they are 
inconsistent with the interpretation of Bitar et al. Thus all numerical results to date are 
consistent with Aoki's proposed symmetry breaking. 

Our second and more important observation is that, close to the continuum limit, one can 
study the pattern of symmetry breaking theoretically, using the chiral Lagrangian. One must 
augment the usual continuum chiral Lagrangian with terms corresponding to the explicit 
breaking of chiral symmetry at non-zero lattice spacing. The form of these terms is dictated 
by the symmetries of the lattice theory. Their coefficients are, however, undetermined aside 
from their order of magnitude. Nevertheless, it turns out that this is sufficient information 
to determine the pattern of symmetry breaking at small but non-zero lattice spacing, up 
to a two-fold ambiguity. If the sign of a particular coefficient is positive, then there is an 
Aoki phase with all the expected properties. If the sign is negative, then there is no Aoki- 
phase, flavor is always unbroken, and the pions do not become massless for any value of the 
bare mass. The present numerical evidence suggests that the first option correctly describes 
(unimproved) Wilson fermions. If this is the case, our calculation makes several predictions 
that can be tested by numerical simulations. 

In Aoki's picture, the flavor-parity violating phases shrink rapidly to isolated points as 
one approaches the continuum limit (see Fig. [I]). At these points, the quarks are massless, 
chiral-symmetry is spontaneously broken, and all pions are massless. Our analysis finds that 
the width of the Aoki phase shrinks as Amo oc a 3 , up to logarithmic corrections. Perhaps 
more importantly, we can see how the flavor-parity breaking at non-zero lattice spacing 
transforms smoothly into chiral symmetry breaking in the continuum limit. In our view, the 
interpretation of Bitar et al. applies only in this limit. 

Similar observations concerning symmetry breaking have been made previously by 
Creutz [[J. He performs a qualitative analysis using the linear sigma-model, and finds the 
same two possible patterns of symmetry breaking at finite lattice spacing. Our calcula- 
tion extends his by providing a firmer theoretical basis for his observations and by making 
quantitative predictions. 

The outline of this paper is as follows. In the next section we review Aoki's proposal 
for the phase diagram, and then describe the alternative view of Bitar et al. In Sec. |T| we 
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explain how the numerical results obtained to date are all consistent with Aoki's proposal. 
We then present, in Sec. fV|, our Chiral Lagrangian analysis. Sec. contains our conclusions. 
We collect technical issues in two appendices. 
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FIG. 1. The phase diagram proposed by Aoki: g is the gauge coupling and mo the dimensionless 
bare quark-mass. The continuum-like phases are labeled A, and the flavor and parity broken phase 
B. The phase diagram is symmetric under mo <-> — (mo + 8). The continuum limit of particular 
interest is that at mo = 0, g = 0. 



We consider QCD with two degenerate quarks, for which the fermionic part of the Eu- 
clidean lattice action is 



where the flavor indices are implicit, and W(mo) is the Wilson-Dirac operator with a (di- 
mensionless) bare mass m common to both flavors. The bare mass is related to the usual 
hopping parameter by mo = 1/(2k) — 4. The action is invariant under 577(2) flavor transfor- 
mations, under parity, and under the interchange mo — > — (m + 8), but it explicitly breaks 
all axial symmetries for any value of mo- 

We are concerned here with the phase diagram in the m -g 2 plane, where g is the gauge 
coupling. Based on perturbation theory, we expect to be able to take a continuum limit by 
simultaneously sending g 2 — > and m — > 0. One can, however, obtain different continuum 
limits by sending mo to any one of the values 0, —2, —4, —6, or —8, for which there are 
respectively 1, 4, 6, 4 and 1 continuum fermions. In the following we will mainly be interested 
in the standard continuum limit at m = (or the equivalent point at mo = —8). If this 
point is approached in such a way that the physical quark- mass is small compared to Aqcd; 
then we expect the theory to spontaneously break chiral symmetry and for there to be a 
degenerate triplet of light pions. This should be true whether we approach from above or 
below mo = 0, since in the continuum the sign of the quark mass is irrelevant. Similar 
phenomena should occur at the other continuum points, except that the symmetry group in 
the continuum limit will be larger (£77(8) at m = —2 and —6, and SU (12) at m = —4), 



II. REVIEW OF PREVIOUS RESULTS 



Cf = ijW(m )ip , 
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and so there should be correspondingly more pions.f] In fact, it may be that the dynamics 
with these large numbers of flavors is different from the tjiq = and —8 cases, since above 
a critical number of flavors we expect there to be no chiral symmetry breaking and no 



confinement [10] 



We are particularly interested in how the phenomena associated with chiral symmetry 
breaking emerge in the continuum limit, given that the lattice theory has no such symmetry. 
In particular, if the continuum limit is approached so that the quarks are massless, then the 
resulting theory should contain massless Goldstone pions. The standard view of how the 
theory accomplishes this is that there is a critical line, m c (g 2 ), along which the triplet of 
pions is massless. This line runs from m = —2 at g 2 = oo to the continuum point at m = 0, 
g 2 = 0. At non-zero g 2 , the pions are not Goldstone particles associated with spontaneous 
chiral symmetry breaking (and in particular their interactions are not restricted by chiral 
symmetry), but they become so in the continuum limit. 



A. Aoki's proposal 



Based on a variety of considerations, Aoki proposed that the vanishing of the pion mass 
is associated with the spontaneous breakdown of flavor symmetry Q. His proposal for the 
phase diagram is sketched in Fig. [I]. The phases labeled A are "continuum-like" in the sense 
that the SU(2) flavor symmetry is unbroken. In phase B, by contrast, flavor symmetry is 
spontaneously broken down to U(l), and parity is also spontaneously broken. There are 
thus two Goldstone bosons in this phase, which are exactly massless in infinite volume, even 
though the lattice spacing is non-zero. We will call these the charged pions. The neutral 
pion which completes the flavor triplet is not a Goldstone particle, and is thus massive within 
phase B. The transitions between the two phases are continuous, so that as one approaches 
the boundary from within phase B, the neutral pion must also become massless so as to 
restore flavor symmetry. This means that as one approaches the boundary from within 
phase A, the degenerate triplet of pions becomes massless. The critical line discussed above 
is the rightmost phase boundary. 

This proposal has been established at strong coupling PJ, for which there are just two 
critical lines (related by the symmetry m — > — (m + 8)). At weak coupling, if one is to 
reproduce the continuum limits expected from perturbation theory, then one must introduce 
the additional regions of phase A shown in the figure^. In particular, if the continuum limit at 
mo = can be approached from positive and negative mass, then the width of phase B must 
vanish at g 2 = 0. Aoki argues that, since phase B is not observed in perturbation theory, 



*At non-zero lattice spacing, however, only an SU(2) subgroup of the full continuum flavor sym- 
metry is exact, so the full pion multiplet will be broken into representations of this subgroup, with 
the lightest pions being in a triplet. 

t These new regions of phase A appear below critical values of g 2 , denoted g 2 x and g 2 2 in the figure. 
These are conventionally drawn as being equal, but we can see no reason to expect this equality. 
The particular ordering of 5^1 an d g 2 2 shown in the figure is an arbitrary choice. 
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its width must be a non-perturbative function of g 2 , i.e. have a power law dependence on 
a. We argue in Sec. [TV] that the dependence is Amo oc a 3 (up to logarithms). 

To investigate the spontaneous symmetry breaking in phase B, Aoki added to the La- 



grangian (2.1) a source term that explicitly breaks both flavor and parity, 

^source = ih^^O^ , (2.2) 

with (T3 acting in flavor space. Flavor and parity are spontaneously broken if, in the limit 
that the source is removed, a condensate remains: 

lim (iil>^5(7 3 ip) = ±const . (2.3) 

Flavor rotations in the 1- and 2-directions change the condensate, while those in the 3- 
direction do not. The symmetry is thus broken from SU(2) — > U(l), with the condensates 
inhabiting the coset space SU(2)/U(1). The condensate ( |2.3| ) appears to violate theorems 
by Vafa and Witten which state that vector-like theories do not spontaneously break flavor 
or parity . As explained in Appendix however, these theorems are not applicable 



to the case at hand. 

Numerical evidence supporting Aoki's proposal has been given in Refs. |l]-|5|. In par- 
ticular, Ref. || studied the pion mass as a function of m at g 2 = 1 ((3 = 6/g 2 = 6) in 
the quenched approximation, finding evidence for all 10 critical lines predicted near the 
continuum limit (although this required the removal by hand of so called "exceptional con- 
figurations"). Unquenched simulations in small volumes also found evidence of Aoki phases, 
although Ref. |3j did not find them at the smallest couplings studied. This could, however, 
be due to the small width of the phases. 



B. The proposal of Bitar, Heller and Narayanan 

The phase diagram for Wilson fermions has recently been reconsidered in Refs. In 
Ref. U, Bitar et al. take issue with Aoki's proposal. They point out that, in the massless 
limit of the continuum theory, a condensate of the form of ( |2.3j ) is simply an axial rotation 
of the usual flavor diagonal condensate, and thus breaks neither flavor nor parity. In other 
words, the condensate ( |2.3| ) spontaneously breaks the chiral group SU(2) L x SU{2) R down 
to an SU (2) subgroup, albeit a different subgroup from the usual flavor group (its generators 
are vector rotations in the 3-direction and axial rotations in the 1- and 2- directions). This 
unbroken subgroup can be defined to be the flavor symmetry. Similarly, the condensate 
preserves a discrete symmetry, which is the original parity transformation conjugated by 
an axial rotation, and this combined symmetry can be taken to be parity. Based on these 
observations, they argue that the Aoki phase should be interpreted as having massless quarks 
and spontaneous chiral symmetry breaking. The condensate of ( j2.3| ) should be interpreted as 
(an axial rotation of) the usual chiral condensate. Flavor and parity symmetry, they claim, 
are not broken at non-zero lattice spacing. 

To study the issue numerically, Refs. |7|,|8| investigate the spectrum of the overlap Hamil- 
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tonian H(—m ) = 75^(^0) -Q The advantage of this operator over the Wilson-Dirac operator 
is that it is Hermitian, and so its eigenvalues are confined to the real axis. This allows one to 
define a spectral density^ p(A; mo) both on the lattice and in the continuum. As explained 
in Refs. the central issue is whether this spectral density, considered as a function of 

A for fixed mo, vanishes for a region around A = 0. In other words, does the spectrum have 
a gap? This is important because, as we now explain, the absence of a gap is indicative of 
chiral symmetry breaking in the continuum. 

We begin by recalling that the Dirac operator in the continuum, unlike the Wilson-Dirac 
operator on the lattice, has a spectrum confined to the imaginary axis. Thus one can define 
a spectral density /5(A) for the operator — ip. As first noted by Banks and Casher [13], chiral 



symmetry breaking occurs if p(A) does not vanish at A = 0. In the continuum, the spectrum 
of the overlap Hamiltonian, which is just 7s(D + m), is related to that of — ip by 



IAI 



VA 2 - m 2 



p(\/A 2 - m 2 ) |A| > 



m 



p(A;m) = <! v - . (2.4) 

|A| < |m| 

We see that at non-zero quark-mass, p(A; m) has a gap between A = — |m| and |m|. As the 
quark mass m — > 0, note that p(A; m) converges to p(A), albeit non-uniformly. Thus the gap 
closes in the chiral limit and p(A; 0) = p(|A|). Based on this continuum result, Ref. f| argues 
that, if the spectrum of 75W / (mo) has no gap, then one should interpret the lattice theory 
as having massless quarks and chiral symmetry breaking. 

As noted in Ref. 0, the spectral density at zero eigenvalue is related to the condensate 
( |2.3|) proposed by Aoki, 

lim (#7503^) = =F2vrp(0; m ) . (2.5) 

h — >0 ± 

Thus Bitar et al. do expect that the Aoki phases should be characterized by a condensate, 
as well as vanishing pion masses. 

References [0,|J study the low-lying spectrum of / ~f 5 W(m Q ) for the mass range —2 < 
mo < 0. They have results on a variety of gauge field ensembles (quenched with and without 
improvement, and partially quenched) all with g 2 m 1. They find that in all cases the gap in 
the spectrum closes for bare masses below a critical line, i.e. for mo < mi(g 2 ). This lies very 
close to the critical curve m c (g 2 ) determined from the vanishing of the extrapolated pion 
mass. We follow Ref. || and assume that the curves would coincide were one to extrapolate 
including the effects of quenched chiral logarithms. What is perhaps unexpected, however, 



^The sign choice of mo is the convention introduced in Refs. [17] and adopted also in Refs. 

§ Spectral densities can be defined on a given gauge configuration by taking the infinite volume 
limit. We refer to these as /?a(A;?71o) in the continuum, and pu(\;mo) on the lattice, with the 
proviso that on the lattice we average over a configuration and its parity conjugate. The quantity 
of interest here, p(A;mo), is obtained by averaging over all configurations, with the appropriate 
weights. For further discussion see Appendix 0. 
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is that the gap does not open up before they reach m = —2. This is inconsistent with the 
results of Ref. JIJ, who find that the Aoki phase has a width of Amo ~ 0.1. 

III. CONSISTENCY 

In this section we address the conflict between the interpretations of Bitar et al. and 
Aoki. We begin with points of agreement. Both contend that there are two phases, and 
that phase B is characterized by a condensate of the form ( |2.3|) . They also agree that the 
masses of all three pions vanish along the boundaries between phases. The conflict concerns 
the number of massless pions within phase B and the interpretation of the condensate. 

We think that Aoki's proposal is a correct description of what happens at any non-zero 
lattice spacing, while the interpretation of Bitar et al. is an approximate description that 
becomes valid in the continuum limit. In this section we elaborate on the former point, while 
in the next we present a framework within which to understand the latter. 

The essential point has already been made by Aoki The only exact symmetry of the 
lattice theory is the SU(2) flavor symmetry, and the condensate ( |2.3| ) spontaneously breaks 
this symmetry down to U(l). Thus Goldstone's theorem implies that there are only two 
massless pions in the broken phase. 

It is worthwhile discussing the symmetry breaking in more detail, so as to emphasize 
the differences between the massless continuum theory and the lattice in phase B. To study 
spontaneous symmetry breaking, one adds a source term that explicitly breaks the symmetry, 
calculates the order parameter (here the condensate), and then takes the limit of vanishing 
source strength. If the symmetry is spontaneously broken from G —>■ H, the resulting 
condensate lives in the coset space G/H. One can thus explore the coset space by examining 
the condensate for sources that effect breaking to different if-subgroups of G. This leads us 
to consider, both on the lattice and in the continuum, the following source: 

£ SO urce(0, h) = hip exp [id ■ <T7 5 ] if) , (3.1) 

where er is a three-vector of Pauli matrices that acts only in flavor space. This source is 
obtained by applying a general chiral transformation to the mass term htpip. We will denote 
the norm of 6 by 9, while will represent a unit vector in the direction of 6. Without loss of 
generality we can take h > and < 9 < tt (so that sin# > 0). In the massless continuum 
theory, standard arguments presented in Appendix [B| lead to 

lim (VvAa) = -7r/5(0) cos 9 5 ab (3.2a) 

h— »0+ 

lim -M)bi>lsipa) = -7rp(0) sm9[0 ■ er] ab . (3.2b) 

h — >0~*~ 

where a and b are flavor indices, and /5(A) = p(A; m = 0) is the spectral density of the 
Dirac operator introduced above. If chiral symmetry is broken down to an SU(2) subgroup, 
then the condensates live in SU(2) L x SU(2) R / SU(2), and should be described by three 
parameters. This is indeed what we find. 

The calculation in the lattice theory is somewhat different, and is explained in Ap- 
pendix 0. Only the pseudoscalar condensate is proportional to the spectral density, and we 
find 
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i- 1 7 ■ ,\ j -7rp(0;m ) [0 • cr] ab 6>0 , , 

ft bm (#75^) = j # = ' ^ ^ 

Equation ( |2.5| ) is a particular case of this general result. The essential difference between 
(|3.3|) and the corresponding continuum result (|3.2b|) is the absence of the factor sin 9. This 



means that the lattice condensate is of fixed magnitude, and varies in direction according to 
the applied source. It depends on only two parameters, the angles in 0, and maps out the 
coset space SU(2)/U(1). This clearly represents a different pattern of symmetry breaking 
than in the massless continuum theory. 

We wish to emphasize that these two different modes of spontaneous symmetry breaking 
may be differentiated by simply enumerating the number of massless pions: in the massless 
continuum theory, there are three massless pions, while on the lattice in phase B there are 
only two. To test this, one needs to know which interpolating fields couple to the different 
pions. For definiteness we consider a source with 6 = n/2 in the 3-direction, i.e. = 
which leads to the condensate {^i^a^ip) ^ 0. This is the choice used in the numerical work 
described in the previous section. Interpolating fields for the Goldstone pions are obtained 
by doing a spatially dependent infinitesimal G/H transformation on the condensate. In 
the massless continuum theory, these transformations are the vector transformations in the 
1- and 2-directions, ip(x) — > exp(—ia(x)ai } 2)^{x) and the axial transformation in the 3- 
direction, if)(x) — > exp(— ia{x)a^)ip{x) . These yield the pion fields 

7Ti(x) = i$(x) (T 2 j5lp(x) , VT 2 (x) = -Z?/)(x)cri7 5 ^(x) , 7T 3 (x) = $ (x)lj) (x) . (3.4) 

On the lattice axial transformations are not a symmetry, and only the first two interpolating 
fields correspond to Goldstone modes. We see that the charged pions are created by the 
familiar pseudoscalar flavor non-singlet operators, although with the flavor indices 1 and 2 
interchanged. The two-point functions of such operators have no disconnected contributions, 
even in the presence of the source term. These are the pions whose masses have been 
calculated in numerical simulations. 

The neutral pion, on the other hand, is created by an operator that appears to be a 
scalar and a flavor singlet. This is because the condensate itself points in a flavor and parity 
violating direction. Indeed, when calculating the neutral pion correlator it is crucial to do so 
with the source turned on, so as to pick out the vacuum state, and then to take the limit of 
vanishing source strength. The correlator for this operator has not been studied numerically. 
We expect such a calculation will be difficult, because it contains both disconnected and 
connected contractions, 

Mx)-K 3 (y))u = lrSu{x; x)JrSu(y; y) - Tr {Su(x; y)Su(y; x)} . (3.5) 

Here Su(x; y) is the fermion propagator in the presence of the source. While in the contin- 
uum, the disconnected contributions can be shown to vanish because of chiral symmetry - 
this must be so since the three pions are degenerate — this is not so at non-zero lattice spac- 
ing. One might hope, however, that their contribution will be small close to the continuum 
limit, and so it would be interesting to calculate the connected contribution alone. In the 
next section we give a prediction for the mass of this particle in the Aoki phase. 
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IV. CHIRAL LAGRANGIAN ANALYSIS 



In this section we study the extent to which the spontaneous breaking of flavor and parity 
at non-zero lattice spacing can be viewed as an approximation to the spontaneous breaking 
of chiral symmetry in the continuum limit. We do so by first determining the effective 
continuum Lagrangian which describes the lattice theory at non-zero lattice spacing, and 
then extracting the effective chiral Lagrangian which describes the long distance behavior of 
this effective continuum theory. The analysis relies solely on symmetries and the assumption 
that dimensionful quantities have a size determined by Aqcd to the appropriate power. In 
the following we shall, for the sake of brevity, drop the subscript on Aqcd- 

Close to the continuum limit, the lattice theory can be described by an effective continuum 
Lagrangian in which the usual terms have been supplemented by contributions proportional 



to powers of the lattice spacing [ H| . These additional pieces are constrained by the sym- 
metries of the lattice action, which in the present case means that they need not respect 
chiral symmetry. The enumeration of operators is identical to that carried out as part of the 
improvement program [^5], and the result is 



+ b2aip(p + m) 2 ip + b3amip(p + m)ip + b^amC g + b 5 am 2 ^ip + 0(a 2 ) , (4.1) 

where C g is the gluon Lagrangian, and m is a physical mass defined below. The symbol 
~ indicates that we are not attempting to control factors of order unity. For example, the 
dimensionless lattice bare mass mo appears in C c g multiplied by a factor Z m (g 2 ,lna) « 1, 
which we ignore. The dimensionless constants m c and bi are functions of g{a) 2 , and in 
general, also of In a. We work at a fixed a and thus suppress this dependence. 

The first two terms in ( [4. 1| ) are the naive continuum limit of the lattice Lagrangian, while 
subsequent contributions result from the fact that the lattice action breaks chiral symmetry. 
The dominant correction at small a is the additive mass renormalization proportional to m c , 
which is linearly divergent as a — > 0. We can combine the two mass terms in the usual way by 
introducing the physical mass m = (mo — m c ) / a (again ignoring Z-factors with logarithmic 
dependence on a). As we will see, m c is very close to, but slightly different from, the critical 
mass m c (g 2 ) at which the pion masses vanish. 

The Lagrangian can be further simplified as follows. By changing quark variables, one 
can remove terms that vanish by the leading order equations of motion [l(J, which in this case 
are those proportional to 62 and 63. Note that the 04-term renders the effective continuum 
coupling constant dependent upon the bare quark-mass, and the 6 5 -term causes the physical 
mass to have a quadratic dependence on the bare quark-mass. While these two contributions 
cannot be removed, they can be ignored. As will become clear shortly, we focus on the 
region in which am ~ (aA) 2 , i.e. physical quark-masses of 0(a). In this region, corrections 
proportional to am are suppressed by aA 1 relative to the terms we keep. With these 
simplifications, the effective Lagrangian becomes 

£ eff ~ £ g + ^(p + m)i) + bia&a^F^ + 0(a 2 ) , (4.2) 

i.e. QCD with a Pauli term. 
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The next step is to write a generalization of the continuum chiral Lagrangian that includes 
the effects of the Pauli term. Without either the mass term or the Pauli term, the theory 
is invariant under SU(2)l x SU(2)r chiral rotations, and its low momentum dynamics is 
described by the chiral Lagrangian 

£2 

£ X = ^-Tr (^E^EJ . (4.3) 

Here E is an SU(2) matrix-valued field that transforms under the chiral group as 

£ -> LEi?t , (4.4) 

with L are R being independent SU(2) rotations. Its vacuum expectation value, So = (E), 
breaks the chiral symmetry down to an SU(2) subgroup. The fluctuations around Eo corre- 
spond to the Goldstone-bosons, 

£ = £ exp | i ^ Tt a a a / f n j . (4.5) 

Adding the mass term to the underlying Lagrangian explicitly breaks chiral symmetry, and 
its contributions to C x can be determined using a standard spurion analysis. Since we are 
only interested in determining the vacuum state, we do not consider derivative interactions. 
Working to second order in m, the potential energy is 



V x = -^Tr(E + E') + £|{Tr(E + Et)} 2 , (4.6) 



with C\ ~ mA 3 and c 2 ~ m 2 A 2 . We have used the property that any term invariant under 
the vector symmetry can be written as a function of Tr (y, + E^ . 

We now include the Pauli term. This is straightforward since it transforms under chiral 
rotations exactly as does the mass term. Thus its effects can be included (along with those of 
the mass term) by making the substitution m — > m + aA 2 . The factors of A are required by 
dimensional analysis, and we recall that we are dropping dimensionless coefficients of order 
unity. Thus the potential is given by fl4.6|), with coefficients 



ci ~ mA 3 + aA 5 , c 2 ~ m 2 A 2 + maA 4 + a 2 A 6 . (4.7) 

Note that choosing am ~ (aA) 2 makes the contributions of the quark mass term and the 
Pauli term comparable. This is as we would expect: if the physical mass is of 0(a), then the 
strength of chiral symmetry breaking due to the mass term should be comparable to that 
due to discretization. 

Since we are keeping terms of 0(a 2 ) in c 2 , we must also consider corrections of 0(a 2 ) to 
the underlying QCD Lagrangian. Fermion bilinears of 0(a 2 ) such as a 2 ipp 3 ip are necessarily 
flavor singlets, and so contribute in the same way as the Pauli term, but suppressed by aA. 
Four fermion operators such as a^ipipipil) can contribute directly to the c 2 term, and give 
an additional contribution proportional to a 2 A 6 . Thus the estimates of ( |4.7| ) remain valid. 
There are of course subleading contributions to these coefficients, suppressed by powers of 
aA. 
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In the following we determine the properties of the chiral theory, and in particular the pat- 
tern of symmetry breaking, as a function of the coefficients c\ and c-i. Since these coefficients 
depend on mo = m c + ma and on a, we can map our results back onto the phase diagram in 
the mo-g(a) 2 plane. We distinguish three regions of quark masses, each successively smaller 
by a factor of aA: 

1. Physical masses, i.e. m/A fixed and small. For such masses, as a — > 0, c\ ~ mA 3 and 
(>z/ci ~ m/A. In this case both the discretization errors and the contribution of the c 2 
term can be ignored, and the symmetry breaking is as in the continuum. 

2. Generic masses of 0(a), i.e. am ~ (aA) 2 . For these, the contribution of the C2 term 
is suppressed by aA, but discretization errors are important in c\. The value of m at 
which the extrapolated pion mass vanishes is shifted by ~ aA 2 (corresponding to a shift 
of 0(a 2 A 2 ) in m ). It is useful to introduce a shifted mass, m! = m — aA 2 , (where, we 
stress, all constants have been set to unity), defined so that C\ vanishes when m! = 0. 
In terms of this new mass the coefficients are 

ci ~ m'A 3 , c 2 ~ m' 2 A 2 + m'aA 4 + a 2 A 6 . (4.8) 

3. Shifted masses of 0(a 2 ), i.e. am' ~ (aA) 3 . For such masses, the coefficients can be 
simplified to 

d ~ m'A 3 , c 2 ~ a 2 A 6 . (4.9) 

The crucial point is that c\ ~ C2, and so the two terms in the potential are comparable. 
As explained below, competition between the two terms can lead to Aoki phases. 

It is immediately apparent from this discussion that the width of the region in which 
new phenomena can occur, due to competition between the two terms in the potential, 
is Amo ~ aAm' ~ (aA) 3 . This is the result announced in the introduction. 

One might be concerned that higher order contributions to the potential will become 
important if the first two terms are comparable to one another. A simple analysis shows 
that the coefficients of the next order terms are generically of size C3 ~ a 3 A 7 . Compared to 
the first two terms, these are suppressed by aA, and can be ignored, except for the possibility 
of a small region of width Amo ~ & 4 i n which all three terms could cancel. Similarly, an 
a 2 A 2 contribution to c\ will shift m! by 0(a 2 ) but will not otherwise change the analysis. 

To determine the condensate S we must minimize the potential energy (|4.6|) . Writing 
X = A + iB • cr with A 2 + B 2 = 1, the potential becomes 

V x = - Cl A + c 2 A 2 , (4.10) 

where the parameter A is constrained to lie between —1 and +1 inclusive. Note that the 
SU(2)v group action, in which L = R in (|4.4| ), leaves A invariant and rotates B by an 
orthogonal transformation. Hence, when the vacuum state So = Aq + iBo • cr develops a 
non-zero B , the flavor symmetry breaks spontaneously to a U(l) subgroup defined by 
exp ji#B • erj, with B being the unit vector in the direction of B . A non-zero value of 
Bq can occur only when \Aq\ is strictly less than one. 
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Based on the discussion above, we treat c<i as a constant of 0(a 2 ), while c\ varies linearly 
with the bare quark-mass. We do not know the sign of c^. It turns out that if C2 is positive, we 
reproduce the properties of the Aoki phase, and so we discuss this case first. The potential 
is then a parabola with its minimum at A m = e, where we have defined the parameter 
e = C\j1ci ~ m'/(a 2 A 3 ). If this minimum lies outside the range —1 to 1, then A is forced 
to one of the boundary points \A \ = 1. This corresponds to S = ±1, and hence the vector 
symmetry SU(2)y is not spontaneously broken. This situation is illustrated in Fig. 0. If 
instead the minimum satisfies \A m \ < 1, then the vacuum is determined by A = A m and 
B ^ 0, as illustrated in Fig. |[ Since B ^ 0, flavor symmetry is spontaneously broken to 
U(l), and the region — 1 < e < 1 thus has the properties of the Aoki phase. 




-10 12 

A A A m 



FIG. 2. The potential V x vs. A for C2 > and \A m \ > 1. The vacuum is at Aq = 1. 




-1 o 1 

A A =A m 



FIG. 3. The potential V x vs. A for c% > and \A m \ < 1. The vacuum is at Aq = A m . 

Let us examine the spontaneous flavor breaking in more detail. In the Aoki phase, the 
direction of B is determined by the source term. To make contact with the numerical 
simulations, we assume a source that aligns B in the 3-direction, and adopt the new pa- 
rameterization S = cos^q + i sin# o"3- The angle #o is determined by e, as explained above, 
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e < -1 

cos6 ={ e -l<e<l . (4.11) 

1 < e 

The vacuum smoothly interpolates between the two flavor-symmetric values. To determine 
the pion masses, we expand £ about the condensate ( |4.5| ), finding 

sin On cos On 9 q , . . 

A = cos 9 - — ^ 7T 3 -— ^ J2< + °«) ■ ( 4 - 12 ) 

« ^^tt a =l 



The potential becomes 



^ (1 - e 2 ) vrf - c 2 e 2 + 0(tt q 3 ) |e| < 1 



r 

and hence the pion masses are 



v*=r- % (N _ 1)E ^ i<h • (413) 



m\ = m 2 2 = o , — = 1 - e 2 for |e| < 1 (4.14a) 



2c 2 
2c 2 



|e| - 1 for |e| > 1 . (4.14b) 



The results are shown in Fig. |J We see that the pions 7ri j2 are the Goldstone bosons of the 
broken flavor symmetry within the Aoki phase, and that all three pions are massless on the 
phase boundaries. 



0.8 - 
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FIG. 4. Pion masses as a function of e for c 2 > 0. The curves are labeled by the flavor of the 
corresponding pion. 

In summary, we have reproduced the phenomenology of the Aoki phases using the chiral 
Lagrangian. In addition, our calculation makes three predictions: 

1. As already mentioned above, for small enough a, the width of the Aoki phase should 
scale as Am ~ a 3 . This should hold up to logarithmic corrections, which we have 
ignored throughout. 
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2. From (|4. 14|) , the mass of the 7r 3 meson within the Aoki phase is predicted in terms 
of the form outside the phase. For example, the slope of m\ versus the dimensionless 
bare mass m should be a factor of two larger as the phase boundary is approached 
from within compared to approaching it from without. In both cases, ~ a when 
e~0(l). 

3. The form of the spectral density of the overlap Hamiltonian at zero eigenvalue can be 
determined. To do this we note that 



4 



ab 



(4.15) 



On the lattice, the l.h.s. is proportional to p(0; mo), while we have found that the r.h.s. 
is proportional to sin^o, and therefore 



p(0;m ) oc VT 



oc m 3 . 



(4.16) 



It should be possible to test some of these predictions. 



We now return to the possibility that the coefficient c 2 is negative. The potential is then 
an inverted parabola, and the extremum A m is an absolute maximum. This means that the 
vacuum state is always at the edge of the allowed range of A, with S = +1 for c\ > and 
S = — 1 for C\ < 0. The potential can be written 

V x = ^ (1 + H) E < + c 2 - |d| + 0(tt o 3 ) . (4.17) 

J tt a=l 

Thus the flavor symmetry is not spontaneously broken for any value of e, and all three pions 
have the same non-zero mass, 




(4.18) 



o 1 — 1 1 1 1 1 — 

-2-10 1 2 

E 

FIG. 5. Pion masses as a function of e for C2 < 0. 
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The previous analysis cannot choose between the c 2 > and c 2 < cases — for this 
we must rely on simulations. As noted in Sec. |T|, all numerical observations to date are 
consistent with there being spontaneous flavor and parity breaking at non-zero lattice spac- 
ing, and hence with c 2 being positive. Thus it appears that an Aoki phase with a width of 
order a 3 will persist until the continuum is reached. There are, however, two caveats that we 
should mention. First, removing the quenched (or partially quenched) approximation could, 
in principle, lead to a change in the sign of c 2 . Second, it is logically possible that, for the 
lattice spacings studied to date, terms of higher order in a are important and conspire with 
the terms we considered above in such a way that the effective value of c 2 changes sign as 
the continuum limit is approached. In this case the Aoki phase would disappear for small 
enough lattice spacing. 

It is interesting to consider in what way the previous analysis changes when one uses 
improved Wilson fermions. The surprising answer is that the analysis is essentially unal- 
tered, even if one removes all errors proportional to a using non-perturbative improvement 
conditions. The reason is that the key discretization error is the 0(a 2 ) contribution to c 2 
[see ( |4.7| )1, and this remains after improvement. Of course, the sign of the coefficient c 2 can 
change, so we cannot predict whether there will be an Aoki phase. What we can say is that 
one or other of the possibilities discussed above will apply.0 This is surprising because in 
both scenarios the non- vanishing pion masses are proportional to a for e ~ 0(1), while one 
would have naively expected the masses to be proportional to a 2 . We see nothing wrong 
with this result, however. In particular, for any fixed physical quark-mass, the discretization 
errors are of 0(a 2 ) for small enough lattice spacing. 

V. CONCLUSIONS 

We have considered the issue of spontaneous flavor and parity symmetry breaking in 
lattice QCD, and its relation to chiral symmetry breaking in the continuum. Contrary to the 
claims of Bitar et al. || , we have shown that flavor-parity breaking can occur at any non-zero 
lattice spacing. Indeed, our chiral Lagrangian analysis shows that such symmetry breaking 
is one of only two options available to the theory as the continuum limit is approached. If 
this option is chosen, as is indicated by present numerical evidence, then our calculation 
shows how the breaking of flavor and parity at non-zero lattice spacing goes over smoothly 
into chiral symmetry breaking in the continuum limit. This is the sense in which we agree 
with the proposal of Bitar et al. that one can think of flavor-parity breaking on the lattice 
as corresponding to chiral symmetry breaking in the continuum. 

Based on this proposal, Bitar et al. suggest obtaining the continuum chiral condensate, 
Co = (uu) + (dd), in the following way. Calculate the condensate (?pij5cr 3 ip) in the infinite 
volume limit at non- zero lattice spacing; average the result over the Aoki-phase; and then take 
the continuum limit. Their major point, with which we agree, is that the lattice condensate 
is a good order parameter, and in particular is free from additive renormalization. It follows 



**In fact the numerical evidence of Ref. Q suggests an Aoki phase also for non-perturbatively 
improved fermions. 
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from our analysis, however, that one should not average over the Aoki phase, but rather take 
the maximum value of the lattice condensate within the Aoki phase, and extrapolate this to 
the continuum limit. If one were to average over the Aoki phase then one would find a value 
for the continuum condensate too small by a factor of 7r/4. 

To understand this, recall that the condensate takes the form E = cos 9 + % sin 9 a 3 in 
the Aoki phase, and so 



Thus at a general position in the Aoki phase, one can obtain continuum condensate from 



On the lattice, however, one can only measure the second of these condensates, and this 
equals Cq only when e = cos#o — 0, i.e. when the lattice condensate takes its maximum 
value. In other words, the lattice condensate picks out a particular component of the chiral 
condensate, and only when the former is maximal are the two condensates aligned. 

One of the predictions of our analysis is that the width of the symmetry-broken phase 
grows rapidly as the lattice spacing is increased, Am ~ a 3 . According to Aoki's picture, 
there is a critical lattice spacing at which the five "fingers" of the Aoki phase merge into a 
single phase of width Am ~ 6 (see Fig. |l|). It is not clear from numerical evidence at what 
coupling this occurs. At g 2 = 1, in the quenched approximation, Ref. |4[] finds evidence for 
all five fingers, while Ref. || finds only a single phase down to mo = —2. We stress, however, 
that either possibility is consistent with our results. Our prediction is that there is some 
critical coupling below which the Aoki phase B splits into finger-like regions of decreasing 
width, surrounded by islands of unbroken phase A, as in Fig. [I]. Once one has entered this 
region then our other predictions, for the pion masses and the condensate, apply as well. 

As explained in Ref. ||, practical applications of domain- wall fermions probably require 
that the Aoki phase has a small width. Our results strongly suggest that this will be the 
case for small enough lattice spacing. 
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APPENDIX A: THE INAPPLICABILITY OF THE VAFA-WITTEN THEOREMS 



The condensate (|2.3|) breaks both flavor and parity, in apparent contradiction with 
Refs. [11] and [12]|, respectively. In this appendix we give a brief explanation of why these 



references do not, in fact, rule out such a condensate. 



The argument of Ref. |TlJ] against flavor-breaking in the continuum does not apply for 
massless fermions because the spectrum of the Dirac operator has a non-vanishing density 
of zero eigenvalues. As has been discussed in Ref. ||, the relevant operator on the lattice is 
the overlap Hamiltonian, 7 5 W / (— m ), and the Vafa-Witten argument also does not apply if 
this has a non-vanishing density of zero eigenvalues, which it does in the Aoki phase. These 
failures are demonstrated by the explicit calculations of Appendix ||. 

The loophole in the argument against spontaneous parity-breaking is simply that the 
analysis says nothing about fermionic bilinears. The bulk of Ref. |12[ is devoted to proving 
that parity-odd bosonic operators do not receive vacuum expectation values. In Footnote 7, 
however, Vafa and Witten claim that the theorem can be extended to include fermionic 
bilinears, since, upon integrating out the fermions, the fermionic operators can be replaced 
by bosonic operators, and for these their original argumentation applies. The example 
they consider is the flavor- singlet operator qi^^q, which is replaced by X = Jri^S a{x; x) , 
with Sa{x',u) being the fermionic propagator in the background A. This extension of their 
argument is not valid, however. The operator X is either ill-defined if there are fermionic 
zero modes, or it simply vanishes because the eigenvalues of 75 (x; y) come in opposite-sign 
pairs. 

Properly speaking, one must calculate the vacuum expectation value by adding a source, 
and only take the source strength to zero at the end. The presence of the source regulates 
the corresponding effective bosonic operator X. However, for two flavors of quarks, it also 
invalidates an essential assumption of Ref. [T2J, namely that X be imaginary in Euclidean 



space. Indeed, one can explicitly calculate the form of X using a spectral decomposition of 
the corresponding fermion operator. For the condensate of interest, ( |2.3| ), for which one adds 
the source ( |2.2|) , this calculation is done in Appendix [FJ. We find that the corresponding 
gluonic operator X is purely real, and so the Vafa-Witten argument does not apply^ This 
point has also been noted in Ref. 0. 

It should be noted, however, that Vafa and Witten's argument does still imply that all 
parity-odd operators composed only of gluon fields do not have vacuum expectation values. 
This suggests that there exists a modification of parity which acts in the same way as parity 



ft For the flavor-singlet parity-odd operator qi^q, on the other hand, the effective bosonic operator 
X in the presence of the source £ S ource = hqexp [iOjs] q is not pure real, and the loophole in applying 
Ref. jl2) is more subtle. Note that the source can be transformed into the mass term hqq by an 
axial U(l) transformation, which, because of the axial anomaly, introduces a OFF term. The source 
therefore does not have the usual effect of picking a particular direction for the condensate, but 
instead selects a different theory characterized by a parity- violating (9-vacuum. This is true even for 
infinitesimal source strength, and thus Ref. [12| does not imply that the fermionic parity-breaking 
condensate (qi^q) vanishes. 
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on the gluon fields, but differently on fermions, and which is unbroken by the condensate in 
the Aoki phase. Such a transformation does exist: it is the product of parity with a discrete 
flavor transformation, 

P' = P®z, z = ia 1 eSU{2). (Al) 

The flavor transformation has no effect on the gluon fields. Note that (P') 2 = — 1 when 
acting on fermion fields, so that P' generates a Z4 group. When acting on the bilinear 
V^75 cr 3' i /'; however, its square is the identity. Note also that other choices of z are possible, 
but all are related by transformations in the unbroken £7(1) flavor subgroup. 



APPENDIX B: CALCULATION OF THE CONDENSATES 

In this appendix we compute the condensates {xp b Yipa) in both the continuum and on the 
infinite-volume lattice with Wilson fermions. We work in two-flavor QCD with a common 
bare mass, and we denote the quark field with flavor index a by ?p a . The quantity T is an 
appropriately chosen spinor-matrix, which in the continuum can be either the unit matrix 1 



or «7 5 , while on the lattice we consider only T = 275. As discussed in Sec. |TJ, the condensates 
are defined by adding to the Lagrangian the source term 

^source h) = hip exp [iO ■ 0-75] ip , (Bl) 

computing {4>bVipa)h f° r non-zero h, and then taking the h — > + limit. We will let 9 denote 
the modulus of 6, while 6 will represent a unit vector in the direction of 6. Without loss of 
generality we can take h > and < 9 < ir (so that sin# > 0). 



1. The continuum theory 

We first consider the massless continuum theory with the fermionic Lagrangian Cf — 
ipapipa. It will be convenient to place the system in a finite box of four-volume V4, so that 
the Euclidean (anti-Hermitian) Dirac-operator ]Z) has discrete eigenvalues ±i\ n with A n > 0, 
and at the end of the calculation the infinite- volume limit followed by the h — >• + limit will 
be taken. For simplicity, let us first work out the flavor diagonal case 6 = z with arbitrary 9. 
After performing the fermionic path integrals, the condensates in the presence of the source 
can be written 

$ b Til> a ) h = JVA exp {-S 9 [A}} [det (JZ> + h)} 2 (^ b T^ a ) A>h , (B2) 

with S g [A] being the pure gluonic Euclidean action, and where the fixed-background con- 
densates take the form 

(ATip a )A,h = -Jr / d 4 x ($ b Tip a (x)) Ajh = --^ Tr T , afc r . (B3) 

V 4 J 14 p + h exp {z6>75cr3 ab \ 

The trace includes a sum over color, spin, and space-time, and it can be written 
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r v 2h i cose for r = i 

r p + /iexp{±^ 75 } ~£z \ 2 n + h? { ± sin^ for T = i-f 5 ' ( > 

where, in computing the fermionic determinant in (|B^ ) and the trace ([B4j) , we have used 
the fact that 



p + m! + i/ii75 has eigenvalues = m x ± iyA 2 + /if (B5a 



75 (p + m 1 + ipi^) has eigenvalues A n = %p\ ± y A 2 + mf . (B5b) 

Since we are working with two- flavor QCD, the fermionic determinant is real, positive, and 
^-independent. 

In infinite volume, the eigenvalues of the Hermitian operator — ip become continuous 
and are described by a spectral density per unit four-volume PaW- Since — ip and 75 anti- 
commute, the eigenstates of — ip come in pairs 4>\{x) and 7s0a(^) with eigenvalues A and 
—A respectively. Therefore, /m(A) is symmetric in A and we can write 

1 2h r°° _ h 
lim lim —/.tt, 7^= lim / c?Apa(A)--^ — r = 7tpa(0) , (B6) 

where we have used the relation 

h-*o+ \±ih A 

with Pj denoting the principal part under the A integral. The principal-part contribution 
vanishes, however, because of the symmetry in A. In the flavor diagonal case we thus find 

lim + (^r^) A , = -7T p A (0) x I ™ 8 *^ T T Z* 15 ■ (B8) 



The condensate in the presence of the source ( |B"T| ) with an arbitrary flavor orien- 
tation 6 can be obtained from ( |B8| ) by an appropriation spinor-rotation. If we write 
= (sin a cos j3, sin a sin (3, cos a) with = (sin(o;/2) cos/?, sin(a;/2) sin/3, cos(a/2)), then 
R = i® ■ (j satisfies 

Ra 3 R j = ■ cr , (B9) 

and we see that R is a spin-1/2 representation of the 3-dimensional rotation from z to 6. As 
75 does not act in the flavor space and since its square is unity, the matrix R also transforms 
the diagonal source with — z to the general source ( |B1| ) using the same transformation as 

i?exp 

[z0<7 3 75] R ] = exp [id ■ 0-75] . (BIO) 

Using (|B10| ), the diagonal condensate flB~8|) can be rotated to a general orientation 6 to 
produce 

(Ma) A = -7r p a(0) cos 6 5 ab (Blla) 

($bil^a) A = -n p A (0) sm9 [6 ■ cr] ab . (Bllb) 
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The condensates (ipbTipa) can now be obtained with the help of ([B2|), and they take the same 
form as (FBI 1|) , except that one uses the gauge averaged spectral density 

p(X) = jvA exp {-S g [A}} [det pf p A {\) . (B12) 

Equation ( |3.2j ) now follows. 

It is also interesting to consider the condensates for a non-zero common mass m with 
fermionic Lagrangian Cf = ip a (p + m) ip a . As with the massless case in (|B6|), a condensate 
before the /i — > + limit has been taken can be expressed as a spectral integral of some rational 
function. The numerators of the scalar and pseudo-scalar condensates are proportional 
to m + h cos 9 and h sin 9 respectively, while the denominators of both condensates are 
A 2 + h 2 + m 2 + 2mh cos 9. The non-zero mass m renders the h —>■ + limit safe for all 
eigenvalues A, thereby producing the ^-independent expressions 

d\p A {X)— -5 ab (B13a) 

-oo A + TYl 

WiMa = • (B13b) 



Of course, taking the m — > limit first would give ( Pll|) , and we see that the m — > and 



h —>■ + limits do not commute. Note, however, that ( |B13a| ) is consistent with the fact that 



flavor is not spontaneously broken for massive quarks [11], since it implies that {^a^) = 



while ( |B13b| ) is consistent with there being no spontaneous flavor or parity breaking. 



2. The infinite-volume lattice theory 

While much of the corresponding calculation on the infinite-volume lattice with Wilson 
fermions is similar to the continuum case, there are some key distinctions which we outline 
in this section. For one thing, the Wilson-Dirac operator W(m ) has complex eigenvalues 
and therefore does not possess a spectral density. This suggests that one tries to relate 
the condensates to the overlap Hamiltonian H(—m ) = ^ 5 W(m ), which is Hermitian and 
therefore has eigenvalues constrained to lie along the real axis. For a fixed background con- 
figuration U, we will write the discrete finite-volume eigenvalues of ■j5W(m ) as A n (m ; U), 
while in infinite volume the corresponding spectral density will be denoted by pj/(A;mo). 
Therefore, the fermionic Lagrangian for the quarks ip a with a common bare mass mo can be 
written Cf = ip a W{mo)%l) a = ip' a H(—mo)ip a , where we now consider ip a and ip' a = ipalh as 
the independent Grassmann variables. While the Wilson term explicitly breaks the chiral 
symmetry, there is nonetheless an exact SU(2) flavor symmetry. We should point out that, 
on the lattice, we only calculate (ijj b Ttp a ) with T = ij 5 . Unlike the continuum, the condensate 
with T = 1 is not proportional to the spectral density, since the eigenvectors of / y 5 W(mo) 
with opposite-sign eigenvalues are not connected by 75. 

We first consider the flavor diagonal case 8 = z with arbitrary 9. As in the continuum, 
we place the system in a finite box of four- volume V4, calculate the condensate with source 
(|B~T|), and then take the infinite- volume limit followed by the h — > + limit. The condensate 
('0&275'0 a )/ l vanishes when 9 = 0, and we therefore concentrate on 9 > 0. After performing 
the fermionic path integrals over if) and ip' = ^75, we find 
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(VWsVUfr = / VU exp{-S g [U] } 



det [7 5 W( m o) + hj 5 exj){i9j 5 }] 



(i>bijsi>a)u,h > ( B14 ) 



with S g [U] being the Euclidean Wilson action, and where the background condensates take 
the form 

$bilMufi = -77- Tr — — — t6ab — r- , (B15) 

V 4 7 5 VV (m ) + hj 5 exp {i9^a 3ab \ 

with the trace including a sum over spin, color, and space-time. Note that in two-flavor QCD, 
the fermionic determinant factors are real and positive, although unlike in the continuum 
there is explicit ^-dependence at non-zero h. Using the chiral basis for the gamma-matrices, 



in which 75 = diag (1,-1), the overlap Hamiltonian takes the general form [17 



H(-mo) = 75^(m ) = ( B + f m ° _ B C _ ^ ) , (B16) 

where the two-component spinor operator C is a discretized chiral Dirac operator, and B 
arises from the Wilson term. From this expression, the operator appearing in the denomi- 
nator of ( |B15| ) can be written 



j 5 W(m ) +/i7 5 exp{z0750- 3o6 } = ^ h W{M) +ihsm6a 3ab . . 

with M = m + h cos 9 , 

and it therefore has eigenvalues A n (M; U) + ihsin9 a 3ab . The infinite-volume limit of the 
trace thus becomes 

J im T7 Tr — W( — VTT, 7^ 1 = / dX Pu(^M) — — . (B18) 

vi^oo V4 7 5 vy (m ) + /i75 exp {i9j 5 a 3ab } J -00 \ + ih sin 9 a 3ab 

Recall that in the continuum, the eigenstates of — ip with opposite-sign eigenvalues were 
connected by 75, and this rendered the spectral density £u(A) symmetric. This is not true 
on the lattice, and the spectral density p[/(A; mo) is not symmetric in A. Nonetheless, if U' 
is the parity conjugate of the configuration U, then p(//(A;m ) = pc/(— A;m ) Hl7| . We can 
thus form a symmetric spectral density by averaging over parity conjugated configurations, 
which we implicitly assume has been done throughout the following. 

All that remains now is to take the h — > + limit. At non-zero lattice spacing, the spec- 
tral density pt/(A;m + hcos9) uniformly converges to p£/(A;m ) as h vanishes, and there- 
fore the limit may be brought inside the integral and allowed to act separately upon each 
term in the integrand. Q As sin^ is strictly positive, it may be absorbed into a redefini- 
tion of h in the denominator of (P18| ), and therefore ( P7|) gives the independent result 



{' l l J b'ilh' l l ) a)u — — t 1 " Pc/(0; mo)cr 3a6 . Remembering that the condensate vanishes when 9 = 0, 
and upon rotating to an arbitrary flavor direction specified by 8, we therefore find 



^In the continuum, there is a problem with uniform convergence for massless quarks, but this will 
dealt with at the end of this section. 
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(falrta)v=[^^™)[0-°U , (B19) 

with the gauge- averaged condensates (V'&rys'^a) taking the same form, except that one uses 
the gauge-averaged spectral density 

p(A;m ) = J VUexp{-S g [U]} [det 75 ^(m )] 2 Ptf(A;m ) ■ (B20) 



Equation ( |3.3|) immediately follows. It can be shown that the gap closes, thereby rendering 
Pu(0;m ) non-zero, only when the dimensionless bare mass m is negative [|T7|]. 

While the calculational techniques presented in this section are not the most natu- 
ral for the continuum, one could nonetheless follow the same formal steps that led up 
to ( P18|) . If the h —>■ + limit can be taken inside the integral, then for quarks of 
physical mass m (and for non-zero 6), one obtains the ^-independent continuum result 
(i J bi'75'4'a) a = — ti"Pa(0; m ) [0 ■ <r)ab, with p^(A; m) being the spectral density for the Her- 
mitian operator 7s(p + m). When m is non-zero, we see from ( |2.4|) that Pa(0; m) = 0, and 
the condensate therefore vanishes. This agrees with with ( |B13b ), which was obtained using 



the methods of the previous section. However, if m vanishes from the start, from (|2.4j ) it is 
apparent that even though p^(A; h) converges to p^(A; 0) as h — > 0, it does not do so uni- 
formly. Thus, we may not interchange the spectral integral and the h — > + limit in ( |B18 ), 



and the calculational techniques of this section break down. In this case, the methods of the 
previous section must be used and (pilb ) is the correct expression for the condensate. 
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